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I. MOTIVATION 



Over the past decade, new definitions of black liole liorizons liave emerged wliicli provide 
powerful tools for studying the behavior of black holes in the strongly dynamical regime. 
These ideas share the common philosophy that black holes should be thought of as physical 
objects in a spacetime that may be identified by local measurements. By contrast, traditional 
black holes and event horizons are globally defined properties of the causal structure of a 
spacetime [1, 2]. 

Though there has always been a certain amount of interest in the dynamics of apparent 
horizons and their relation to black hole physics (see for example [3]), Hay ward began this 
line of research in earnest with his definition of trapping horizons. These were initially 
used to formulate dynamic versions of the laws of black hole mechanics in a quasi-local 
context [4]. Since then however, they have found a variety of applications including, for 
example, studying interactions between black holes and gravitational waves [5, 6]. A bit 
later, the isolated horizons of Ashtekar et al. [7, 8] were developed to provide a quasi-local 
characterization of the equilibrium states of black holes. In those works, it was shown that 
these objects obey a phase space formulation of the zeroth and first laws of black hole 
mechanics. Further, loop quantum gravity has made important use of them as boundary 
conditions in calculations of black hole entropy [9] . 

Closely related to both trapping and isolated horizons are dynamical horizons [10, 11], 
which characterize the dynamical phase of smooth black hole evolutions. In particular, it has 
been shown that flux laws can be formulated for these horizons that measure the growth of 
such quantities as energy and entropy. These laws contain terms that may be identified with 
fluxes of particular physical quantities such as matter and gravitational waves. A similar law 
exists for trapping horizons [12] and it has been shown that the energy expressions involved 
agree with those derived in a Hamiltonian analysis of horizons as spacetime boundaries [13] . 

In other developments, studies have been made of the perturbative, "almost isolated" 
regime [5, 14, 15] and isolated horizons have been given a convenient characterization in 
terms of multipole moments, which can be extended to a definition for the multipole moments 
for dynamical and trapping horizons [16]. Mathematical investigations have been made into 
such properties as existence and uniqueness [17, 18] and recently these notions have also 
begun to find application in the physical interpretation of numerical results [19-21]. 

Thus, isolated, dynamical and trapping horizons constitute an increasingly well-developed 
quasi- local framework for analytical as well as numerical studies of black hole dynamics in the 
strong field regime. However, in order to make effective use of them it is important to have 
clear intuition about how they behave. Unfortunately, there is a paucity of good, analytical 
examples of spacetimes containing dynamical and/or trapping horizons and this has recently 
given rise to some confusion about the generic behavior of dynamical and trapping horizons. 
In this paper we will attempt to clarify matters by presenting several analytic and numerical 
examples of spacetimes containing these horizons. 

Before considering these in more detail, let us recall the relevant definitions. From 
Hay ward [4] we have: 

Definition 1. A trapping horizon H is a hypersurface in a 4-dimensional spacetime 
that is foliated by 2-surfaces (which we will take to be of spherical topology) such that 
d{i)\H = 0, 9(ri)\H 7^ and CnO(£)\H 7^ 0. A trapping horizon is called outer if Cn0{£)\H < 0, 
inner if £„6'(£)|// > 0, future if 9(n)\H < and past if 0(^n)\H > 0. 



2 



In this definition, and what follows, £° and are respectively the future-directed 
outgoing and ingoing null normals to a leaf of the foliation while and are the 
expansion of the congruences of curves generated by those vector fields. Further, it is 
assumed that has been extended so that it is surface generating (ie. n A dn = 0) in 
some neighbourhood of H. Then, from the definition, outer trapping horizons have trapped 
surfaces "just inside" them while inner trapping horizons have trapped surfaces "just 
outside" . 

We will mainly be interested in a slight generalization of future trapping horizons that 
was recently introduced by Ashtekar and Galloway [17]: 

Definition 2. A marginally trapped tube (MTT), T is a hypersurface in a 4-dimensional 
spacetime that is foliated by two-surfaces (again assumed to be of spherical topology) such 
that 0(^n)\T < and 6'(^)|t = 0. 

We refer to the leaves of the foliation as marginally trapped surfaces} MTTs have 
no restriction on their signature, which is allowed to vary over the hypersurface. However, 
if an MTT is everywhere spacelike it is referred to as a dynamical horizon [8, 11], if it 
is everywhere timelike then it called a timelike membrane (TLM) [11, 17], and if it is 
everywhere null and non-expanding then we have an isolated horizon."^ 

Note that the distinction between dynamical horizons and timelike membranes is more 
than just a technical difference of signature. For example, it is clear that since dynamical 
horizons are spacelike they may only be crossed in one direction by causal curves; this is a 
key characteristic of black hole horizons. By contrast, a timelike membrane, being timelike, 
obviously does not share this property. Further, it may be shown that while dynamical hori- 
zons always expand as they evolve, timelike membranes always shrink. Finally, in contrast 
to the dynamical horizon flux laws of [8, 11], both the geometrical and matter contributions 
have indefinite signatures in the corresponding equations on timelike membranes. In par- 
ticular, the geometric term no longer has a natural interpretation as a flux of gravitational 
radiation energy. 

Simple explicit examples of trapping horizons/MTTs are provided by the Vaidya space- 
times [11]. These are spherically symmetric solutions to Einstein's equations describing the 
formation/growth of a black hole as null dust falls in from infinity. In the absence of a 
cosmo logical constant, the growth of the resulting black hole is always characterized by a 
dynamical/future outer trapping horizon. In the presence of a cosmo logical constant the 
situation is a little more complicated as a second MTT appears — a timelike membrane 
which can be associated with the cosmological horizon. Still, even in this case, the MTT 
associated with the black hole remains a dynamical horizon. 

Heuristic arguments presented in [11] and [22] suggest that under physically reasonable 
circumstances the MTTs associated with black hole formation and growth will always be 



^ This follows the terminology used in [17]. Note however that the exact definition of marginally trapped 
surfaces varies somewhat in the literature. For example, Wald [2] defines a marginally trapped surface to 
be a compact, spacelike, two-surface for which both null-expansions are non-positive. 

^ More precisely, this is a non-expanding horizon of which isolated horizons are a special case. However, we 
will follow the common usage of isolated horizon in its general rather than specific meaning throughout 
this paper. 



3 




FIG. 1: A schematic of an MTT as informed by our later examples. In this diagram 45° lines are 
null, time increases in the vertical direction, while horizon area increases as one moves to the right. 

spacelike and so future outer trapping horizons are generic in this context. This intuition 
needs to be amended. In Oppenheimer-Snyder spacetimes [25], which describe the gravita- 
tional collapse of spherical dust clouds, it has been shown that timelike membranes, rather 
than dynamical horizons, appear during the formation of black holes. These spacetimes are 
constructed out of a piece of the closed FRW universe (which forms a homogeneous and 
isotropic dust ball) surgically inserted into a Schwarzschild spacetime. As this "star" is 
allowed to collapse the horizon structure develops in the following way. At first there are no 
horizons. Then, as the dust reaches a critical density an isolated horizon forms that cuts off 
the continuing collapse from the rest of the universe (in this case it is also an event horizon). 
Coincident ally, a timelike membrane appears and contracts inside the isolated horizon until 
it reaches zero area. 

In this case, the timelike membrane is clearly associated with the formation of the black 
hole and cannot be dismissed as "cosmological" . One might argue that the OS spacetime is 
not very physical. Nevertheless, the question remains as to what is the generic behavior of 
MTTs during the formation or growth of a black hole. When are they spacelike and when 
are they timelike? Equivalently, when do they grow in the way that one would intuitively 
expect and when do they shrink? 

In this paper, we present several new examples of MTT-containing spacetimes. From 
these, we see that in most circumstances horizons are either isolated or dynamical, spacelike, 
and expanding. However, under certain conditions more interesting behaviours are seen. 
There are two ways in which new horizons may form where they did not exist before. In 
the first, an MTT may appear out of a singularity (such as a point of infinite density). In 
the second, large (though not infinite) concentrations of matter may force the creation of 
a dynamical horizon-timelike membrane pair. Such formations are closely related to the 
well-known phenomena of "jumping" apparent horizons [1]. 

Horizons may also disappear. Timelike membrane-dynamical horizon pairs can annihilate 
each other or alternatively MTTs can vanish into singularities. Figure 1 schematically 
displays some of these behaviours, and also suggests another interpretation. One can think 
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of a single MTT that winds its way forwards and backwards through time rather than 
multiple dynamical horizons, timelike membranes, and isolated horizons that appear and 
disappear. 

For simplicity we will usually only consider spherically symmetric spacetimes. Manifolds 
will be smooth unless stated otherwise. We use units such that G = c = 1, and spacetimes 
have signature (— , +, +, +). We will assume that the Einstein equations hold, with matter 
satisfying the null energy condition. 

The paper is structured as follows. In section II we explain the basic properties of 
MTTs. In particular, we discuss how the signature of the metric induced on an MTT by the 
spacetime metric depends on the matter present and then examine several different types of 
matter. In section III we illustrate the analytic evolution of MTTs in various situations for 
the case of pressureless dust. In section IV we consider a couple of numerical examples of 
MTTs in the presence of a massive scalar field. Conclusions are presented in section V. 



II. MARGINALLY TRAPPED TUBES 
A. General properties 

As we have seen, the definition of MTTs is weaker than that of both dynamical horizons 
and timelike membranes (there is no restriction on the signature of the metric on T) , and of 
trapping horizons (there is no restriction that Crfi{t) 7^ 0). On an MTT, both the induced 
metric signature and Crfiii) are allowed to vary. 

There are simple procedures which may be used to determine this signature. Here we 
will restrict ourselves to the case of spherically symmetric spacetimes; the general case is 
closely related (see, for example, [4] or [11]). With this condition, it is natural to restrict 
our attention to similarly symmetric MTTs. ^ Further, it will be convenient to use the 
symmetry to foliate the spacetime into spacelike two-spheres and so extend the definitions 
of and Qi^t) off of T. We require that all such quantities share the symmetry and for 

simplicity will also impose the standard condition that I ■ n = —\. 

Next, following the conventions of [13, 15], we let V" be a vector field which is: 1) 
tangential to T, 2) everywhere orthogonal to the foliation by marginally trapped surfaces, 
and 3) generates a flow which preserves the foliation. Thus, if w is a foliation label, L-\^v is 
a function of v only — it is independent of the exact position on a leaf. Then it is always 
possible to find a function C and normalization of 1°" such that V" = — Cn°- . Moreover, 
the definition of V"' implies that C\iQ{^t) = 0, which gives us an expression for C: 

C = (2.1) 

Note that VaV" = 2C, so that the sign of C determines the signature of T: if C > it is 
spacelike, if C = {CgO^tj = 0) or becomes undefined {CeO{e) 7^ while CnO{e) = 0) it is null, 
and if C < it is timelike. In addition, the sign of C determines whether T is expanding, 



^ Non-spherically symmetric MTTs will also exist in these spacetimes. However, they are more complicated 
to locate and so will not be considered in this first analysis. See [31] and references therein for a discussion 
of what is know about the allowed behaviours of such MTTs. 
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contracting, or unchanging in area. Indeed, denoting the area element of the two-sphere 
cross-sections by e, one has 

Cvi=-Cd(^n)~e- (2.2) 

This means that the expansion or contraction of an MTT is hnked to its signature (since 
9{n) < 0). In particular, when T is spacelike it expands while when it is timelike it contracts. 

Further (2.1) explains the close relationship between the various trapping horizons and 
MTTs. If the null energy condition holds, then the numerator is non-positive (by the 
Raychaudhuri equation). Thus, the sign of C is determined by the sign of CnQ(i)- Away 
from isolation and cases where CnQ(i) = 0,^ an MTT is a dynamical horizon if and only if 
it is a future outer trapping horizon. Similarly it is a timelike membrane if and only it is a 
future inner trapping horizon. 

To understand the generic behavior of these spherically symmetric MTTs we focus on 
this function C. Then, keeping in mind that 6(^e) = 0, from the Raychaudhuri equation it 
is easy to see that Ce9(e) = — 47rTab£"£'' while from Gab^""n^ = SnTabi'^n'' one can show that 
(^nG(tj = —TZ / A + ATiTah(-""n}' , where IZ is the scalar curvature of the two-sphere cross-section 
of the MTT foliation (the easiest way to see this is to consult a table of the Newman-Penrose 
equations, for example in [26, 27]). Then, using eq. (2.1), 

where A is the area of a two-sphere cross-section of the MTT. As noted above, if the null 
energy condition holds then the numerator is non-negative. Thus, the sign of C depends on 
the relative magnitude of 1/(2^4) and Tabi°'n''. 

The MTT behaviours implied by (2.1) and (2.3) are closely related to Theorem 2 of [18], 
which may be summarized as follows. Suppose a (not necessarily spherically symmetric) 
spacetime is foliated by a family of spacelike hypersurfaces. Let S* C Sq be a marginally 
outer trapped surface, i.e., 6'(£) = for an outgoing null normal while the ingoing null 
expansion is unrestricted. Suppose S is also "strictly stably outermost", which roughly 
means that if 5* is deformed outward, the corresponding deformation of the outgoing null 
expansion is non-negative and positive somewhere. Then first of all, S is contained in a 
horizon H foliated by marginally outer trapped leaves that lie in S^, which exists at least as 
long as these leaves remain strictly stably outermost. Moreover, if the null energy condition 
holds, H is achronal. If Gab^"'^^ > somewhere on S, then H is spacelike everywhere near 
So. 

The condition that S be strictly stably outermost is equivalent to the requirement that a 
certain operator Ls„ acting on functions ip on S has a strictly positive principal eigenvalue. 
Restricting ourselves to spherical symmetry again and using the Einstein equations, this 
operator reduces to 

Ls^V- = -A^ + (^^n - 87rT,,£V^ t/^, (2.4) 

where A is the Laplacian operator on the round 2-sphere S and TZ = Sir /A the scalar 
curvature of S. The principal eigenvalue (corresponding to ip = constant) is then 

X = 87r[l/{2A)-Tabrn']. (2.5) 



To understand why the latter imphes an important distinction, see the example in [24]. 
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According to the theorem in [18], if the null energy condition holds then A > implies 
local achronality of the horizon, consistent with our considerations above. If, moreover, 
GabC-""^^ = SnTabi"'^'' > on S" then the horizon must be spacelike, which is again as we 
found. Thus, the main results of [18] applied to the case of spherical symmetry are neatly 
encapsulated in the expression (2.3) for the function C, which, however, will also tell us 
when we are dealing with a timelike membrane. 

We now consider eq. (2.3) for some particular types of matter. 

B. Behavior for some matter sources 

1. Timelike perfect fluid 

For a perfect fluid that moves along timelike worldlines with unit tangent m", the stress- 
energy tensor takes the form 

Tab = (P + P)UaUb + Pgab , (2.6) 

where p is the matter density of the fluid and P is the pressure. Writing = + {2^)~^n°' 
for some function ^, we find that 

^_ 1 P+P 

2e {i/A) + p-p' 

and see that a priori, the MTT could have spacelike, null, or timelike evolution, the deciding 
factor being the magnitude of (p — P) relative to 1/A. 

Simple examples of the various behaviours may be found in Robertson- Walker spacetimes 
[23, 28]. In the case where these cosmological models are collapsing, one can find spherical 
MTTs through all points in M. Picking a single MTT for definiteness (or equivalently 
selecting a "centre" for the universe), and assuming an equation of state of the form P = ap, 
where a is some constant, these surfaces are: i) timelike and contracting if a < 1/3 (this 
includes a dust-filled universe, a = 0, and so the timelike contractions seen in Oppenheimer- 
Snyder collapse), ii) null and contracting if a = 1/3 (a radiation filled universe with divergent 
C), and iii) spacelike and expanding if a > 1/3.^ For now, we are not too concerned with 
the physical interpretation of such models and evolutions, but instead just note that these 
are all allowed mathematical possibilities. 

In the case of pressureless dust, P = 0, one has the following heuristic picture deciding 
how the MTT will behave at a given two-sphere cross-section. Foliate spacetime by the level 
surfaces of dust particle proper time. Fix such a slice, let r be some radial coordinate, and 
define the "mass enclosed by a sphere with coordinate radius r" to be m{r) = moN{r), with 
rrio the mass of a single particle and N{r) the number of particles with radial coordinate 
smaller than r. Then eq. (2.7) can be rewritten as 

C = (2.8) 

PT- P 



^ Of course, these potential behaviours are not confined to barotropic equations of state. In general, given 
the dominant energy condition and assuming that p > P, the signature of the MTT is determined by 
p-'iP. 
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(2.7) 



where / = i?p/(4^^) > with R = {A/AnY^'^ the areal radius, and pt and p are two 
quantities with dimensions of surface density, p = dm/dA is the change of enclosed mass 
with area, while px is a purely geometric quantity associated with the MTT: 

Pt = ^MtTZ, (2.9) 

with Mt = R/2 the instantaneous mass^ of the MTT and TZ = Sir /A the scalar curvature 
of the 2-sphere cross-section. pT is a straightforward generalization of the mass aspect for 
isolated and dynamical horizons [16]: it determines the mass multipoles of the MTT and as 
such its definition and physical meaning are not restricted to spherical symmetry. Thus, the 
relative magnitude of the "matter surface density" p and the "geometric surface density" 
Pt is what governs the behavior of the MTT. Given a 2-sphere cross-section of the MTT, if 
the matter surface density is larger than the geometric surface density then the MTT will 
contract because of eq. (2.2), and in that case it must be timelike. If on the other hand the 
matter surface density is smaller than the geometric one, the MTT will expand, in which 
case it must be spacelike. When the two balance each other the MTT is null (though not 
isolated, as in this case C ^ oo rather than 0). 



2. Null fluids 

We next consider a null fluid that moves inwards from infinity towards some centre with 
tangent vector n". Then, the stress-energy tensor is 

Tab = {p + PWrif, , (2.10) 

and so one quickly finds that 

C = 2A{p + P)>0. (2.11) 

Hence, for matter of this type, the MTT can only be either isolated (if p+P = 0) or spacelike 
and expanding (otherwise). Indeed, this is the type of matter in the Vaidya spacetime where 
it is well known that MTTs demonstrate only these behaviours [11]. 

We note in passing that the Vaidya spacetimes can be generalized to include both outgoing 
and ingoing null dust, plus a distribution of energy whose rest frame is stationary. In such 
cases, the MTT can again be timelike, as follows from the discussion in [29] . We suspect that 
with a careful choice of the parameters in these models, one could also generate examples 
of MTTs that are partly timelike, partly null, and partly spacelike. However, we will not 
discuss these examples further here, choosing instead to focus our attention on the more 
astrophysically relevant timelike dust spacetimes of section HI. 



3. Scalar fields 

The final matter model that we will consider is that of a scalar field 0. This has stress- 
energy tensor 

Tab = ^ ha<P^b(p - ^V,0V^0^7a, - V{(P) gab^ , (2.12) 



^ See [11] for a detailed discussion in the context of dynamical horizons; similar considerations hold for 
MTTs. 
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where V{(j)) is a potential function. For arbitrary potentials, C takes the form 



C - (213) 

If = (a massless scalar field), then £^0 ^ implies that the MTT must be spacelike. 
However, unless V is negative definite, a non-zero potential a priori allows for spacelike, 
null, and timelike evolutions. In particular, this is the case for a massive Klein-Gordon field 
where V{(j)) = mo0^/2 for some rrio > 0. 



III. PRESSURELESS DUST 

Examples of all of the MTT behaviours discussed above can be seen within the Tolman- 
Bondi family of solutions. In this section we will review this surprisingly rich set of solutions 
and then discuss specific members which display the various behaviours. 



A. Tolman-Bondi solutions 

These solutions describe the gravitational collapse of spherically symmetric dust clouds. 
They are very easy to work with and allow us to trace the evolution of a spacetime from 
specified initial conditions. A nice discussion can be found in the early part of [30] and with 
minor changes, we will follow that description below. 

Initial conditions are given on a spherically symmetric, spacelike, three-surface Sq. On 
that surface, we may specify: i) the dust density po, so that on that surface Tab = PoUaUb 
where u"" is the forward-in-time pointing timelike unit normal to S^, and ii) the initial areal 
velocity Vo = ^ of the dust, where r = A/Att is the areal radius and r is the proper time 
as measured by observers comoving with the dust. Then, taking the areal radius r as a 
coordinate on E along with the usual spherical coordinates {0,(j)), these two functions are 
sufficient to specify both the intrinsic metric hab and extrinsic curvature Kab of Sq. Defining 

m{r) = An po{f)f^df and (3.1) 
Jo 

kir) = ^-^-vlir), (3.2) 

we have 

ds^= ' , + r\de^ + smHd<P^) , (3.3) 
1 — k{r) 

and 

Kab = -1-rarb H , (3.4) 

dr r 

where Qab = [d0]a[d9]b+sm'^ 6[d(l)]a[d(f)]b and is the spacelike, unit normal, outward pointing 
radial vector. With the restriction that So be spacelike, we note that initial conditions must 
be chosen so that k{r) < 1. 

The functions m{r) and k{r) have well-defined physical interpretations: m(r) is a mass 
function which measures the amount of matter contained within a sphere of areal radius r 
while k{r) determines whether or not the system is gravitationally bound. We restrict our 
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attention to gravitationally bound systems {k{r) > 0) and so the allowed values of k{r) are 
< k{r) < 1. Note too that if Vo{r) = 0, So will represent an instant of time symmetry 
with Kab = 0. Then k{r) = ^^^^ and so < 2m(r) < r. 

If we further restrict our attention to initial conditions for which all matter is initially 
either stationary or infalling {vo{r) < 0), it is not hard to use the Einstein equations to show 
that with r as the proper time measured by observers comoving with the dust, 



,,,,.^._^|iirl_.,, ,3.5, 

where R{t, r) is defined as the areal radius at time r of the dust shell that had initial areal 
radius r on Sq. 

Then, in Gaussian normal coordinates, these initial conditions can be evolved to give us 
a four-dimensional metric: 

ds^ = -dr^ + + R\r, r) dn' , (3.6) 

1 — k{r) 

where R' = dR/dr and the stress-energy tensor takes the form Tab = p{T,r)VaT^bT with: 

( ) = 1 ^ repair) , . 

^^^'"'^ 47iR^{T,r) dR R^{T,r)R'{T,r)' 

Now, there is an exact, parametric, solution to equation (3.5). Specifically, for < r/ < tt 
and initial areal radius r: 

m(r) 

r(?7,r) = 'ro[r) + j^y^-^{v + smr]) and (3.8) 

^(-)-^-^(i)' 

where 



, , rvn(r) m(r) / rv'^(r) 

rJr = — ^ + — ^arccosWl 3.10 

^ ^ 2k{r) k^l^ir) \j 2m(r) ^ ' 

In the special case where Vo{r) = 0, then To{r) = and r = corresponds to 77 = 0. For 
simplicity, our explicit examples will be restricted to such evolutions; it turns out that these 
are sufficient to demonstrate all of the potential MTT evolutions. The equations (3.8) and 
(3.9) then reduce to: 



R{ri,r) = rcos^ . (3.12) 

In constructing our examples, we will sometimes find it convenient to excise the interior 
or exterior part of a dust spacetime and replace it by a Schwarzschild region. Then, the 
Einstein equations are satisfied at the three-dimensional junction surface if and only if its 
intrinsic and extrinsic curvatures are the same whether measured on the interior or the 
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exterior of the surface. For our purposes, it will be sufficient to only consider excisions along 
the comoving surfaces of constant r. Specifically, suppose we wish to take the junction 
surface as r = f, making the spacetime Schwarzschild either for r < f or r > f. Then 
assuming Vo = 0, the induced metric from the dust side on the corresponding 3-dimensional 
timelike junction surface is 



■^3 

^l + cosr/)2dr/2 + f2cos^ (^1) (rfe2 + sin2e#2)_ ^g^^g^ 



6m(r] 



Kee = ^ = -f cos^ ( ^ ) ^/ 1 - (3.14) 



The induced metric from the Schwarzschild side is the same with m(f) replaced by the 
Schwarzschild mass M . The non-zero components of the extrinsic curvature from the dust 
side are 

^ = -"°^ (2 

and from the Schwarzschild side we again get the same with m{r) replaced by M. Thus 
the resulting restrictions are quite simple. If it is the exterior that is being replaced by 
Schwarzschild beyond some coordinate radius f, then the matching of intrinsic and extrinsic 
curvatures of the boundary requires that the Schwarzschild mass of the exterior geometry 
be M = m(f). If the interior is being replaced by a Schwarzschild geometry with mass M 
for r < f, then the relationship between mass function and initial density should be 

m{r)=M + ATT f^po{f)df (3.15) 
J f 

for r > r. 

Thus, given initial conditions, we can analytically calculate the full, four- dimensional 
metric as long as the Gaussian normal coordinate system remains valid. As can be inferred 
from the preceding discussion, the reason for this simplicity is that any shell of radius r 
effectively evolves along the geodesies of a Schwarzschild solution of mass m{r) — outer 
shells do not affect the evolution of inner shells and inner shells only affect outer shells 
through the total mass function m(r). This makes these spacetimes especially easy to work 
with but it also points to a potential problem. If the shells of constant r do not maintain their 
original ordering, then the mass contained within a shell can change with time, and in this 
case the evolution will no longer be described by the solution discussed above. Further, apart 
from this physical problem, shell-crossings will also cause the Gaussian normal coordinates 
to break down as those surfaces of constant r were also used as coordinates. 

There is a large body of literature on shell crossings (see for example the references in [30]), 
and the Tolman-Bondi spacetimes are often used to study this phenomenon. In this paper, 
however, we are only interested in using these spacetimes to provide concrete examples of 
MTT evolutions. As such we will studiously avoid such complications. A sufficient condition 
to guarantee that no shell-crossings occur is easily seen to be 

i?'(r,r)>0, (3.16) 

as this will ensure a physical separation of r = constant surfaces. If Vo{r) = it is not hard 
to see that this reduces to 

dTc 

->0, (3.17) 
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where 




^c(r) = 7r^-^, (3.18) 
8m r 



is the time for the shell of initial areal radius r to collapse to zero area. 



B. MTTs in Tolman-Bondi spacetimes 

The above considerations define the Tolman-Bondi spacetimes. We now consider the 
location of spherically symmetric marginally trapped surfaces within those spacetimes. For 
definiteness we choose 

R'(t r) 1 R'(t r) 

L = -[drl + ^^^[drl andn, = --[drU / ' ; Idr^, (3.19) 

— k{r) ^ 2a/1 — k[r) 



as the forward- in-time outward and inward null vectors satisfying i ■ n = —1. Note too 
that these expressions will always be well-defined since our requirement that Eo be purely 
spacelike fixes k{r) < 1. 

Of course £ and n can be rescaled, but for the purposes of this paper that scaling is 
irrelevant (we only care about whether quantities defined with respect to them are zero, 
negative or positive) and so we can choose to work with this convenient choice of vectors. 
Then, keeping in mind that 6'(£) = {g"-^ + £°-n^ + n°'t)Va£b (and similarly for 9{n)): 



2(%,r) + v/r3fcH) . n %,r)-v/r^%) 
= ^(.) = ^ . (3.20) 

Solving = with a bit of help from evolution equation (3.5), it is not hard to see that 
marginally trapped surfaces occur whenever 

i?(r,r) = 2m(r). (3.21) 

This is the expected result if one recalls that shells of constant r essentially move along the 
geodesies of a Schwarzschild spacetime with mass m(r). Furthermore, it is clear that on 
such a surface, 

Oin) = -- rT^<0. 3.22 

m[r) 

Thus, in these spacetimes, any two-sphere on which 6'(£) = is part of an MTT. 

Finally we can calculate the expansion parameter C. Using condition (3.21), definition 
(2.1), the evolution equation (3.5), and definition (3.1) one can directly calculate 

^ ^ 2m'{r) ^ 2p(r,r) 

R'{T,r) -m'{r) l/A-p(r,r)' ^' ' 

where A = 167rm^(r). As would be expected this result agrees with the earlier, more general, 
perfect fluid results (2.7) and (2.8). 

With this background established, generating examples of MTT spacetimes is simply 
a matter of picking initial conditions and using (3.8) and (3.9) to generate the evolution. 
The expression for C, eq. (3.23) can then be used to determine the signature of the tube. 



12 




c) T vs R for MTT evolution 
FIG. 2: Collapse of a dust ball 

Following these simple procedures we generate the examples below, all of which will start 
from an instant of time symmetry [vo = 0). Often we will consider situations where the dust 
accretes onto a pre-existing hole. In those cases, excisions will be performed inside of some 
shell so that a black hole may be inserted into the spacetime. When this is done we will use 
the Schwarzschild mass of the interior as a reference scale for masses, lengths, and times. In 
other cases the scale will be set according to the physics of the particular situation. 



C. Dust ball collapse 

1. Collapse of a dust ball with Gaussian initial density 

We begin with the example of a non-uniform dust ball collapsing to form a black hole. 
The initial density distribution is taken to be radially Gaussian so that 

P'^^-^-'"'- (3-24) 
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r = 




T = r] = (reflect through this axis) 



FIG. 3: Penrose- Carter diagram for dust-ball collapse. To obtain the full spacetime diagram, 
including the corresponding white hole and reflect this diagram through the r = r/ = 

instant of time symmetry. Note that r = starts as a regular timelike geodesic in the centre of 
the ball and then changes into a spacelike singularity at r = Tc(0) as the central density goes to 
infinity. 

In this expression rrio is the total mass of the cloud as read from the corresponding mass 
function m(r) (equation 3.1) and also corresponds to its ADM mass. The parameter ro 
determines how much the cloud is initially "spread out". Such a distribution is pictured in 
figure 2a), where we have taken = 100 p is plotted in units of and r in units of 
nio- 

For this choice of parameters r — 2m(r) > everywhere and so there are no marginally 
trapped surfaces in the initial time slice. Further, for any distribution of this kind condition 
(3.17) is met and therefore there are no shell crossings. In this case. 



and the first r = constant shells collapse to zero area at that time. As this happens, the 
central density p(rc, 0) diverges to infinity (equation 3.7). Referring to 2c) we see that an 
MTT is born out of this divergence, while 2b) shows that it is everywhere spacelike and so 
is a dynamical horizon. 

Note that in figure 2c) (as well as in future evolution graphs), the MTT is shown as 
the thick black line. The rest of the lines are r = constant surfaces and correspond to the 
timelike geodesies that trace the evolution of individual dust shells. In this figure, most of 
them are very nearly horizontal as they have fallen in from a long way out and are moving 
very quickly (relative to the constant r foliation) by the time that they approach the horizon. 

Returning to the evolution graph 2c) we note that as r ^ oo, -R — 2mo and C — as 
the last bits of matter fall through the horizon and it asymptotes towards a null and isolated 
state. 

Many of these features of the evolution may also be seen in the Penrose-Carter diagram 
for the spacetime which is given in figure 3. Thus, we again see the MTT created out of the 




(3.25) 
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c) T vs R for MTT evolution 
FIG. 4: Approach to Oppenheimer-Snyder collapse 

central spacelike singularity, evolving in a spacelike fashion, and finishing its evolution by 
asymptoting to the null event horizon. Note that for simplicity of presentation, this diagram 
shows a sharp cut-off of the dust distribution at some finite r. In our example this does 
not occur and instead the density asymptotes to zero. Thus, properly one should view the 
cut-off as as marking, say, the r which contains 99.9% of the mass. This diagram also nicely 
shows that while the MTT appears at the same time as the singularity, nothing in particular 
is happening at r = when the event horizon "appears" . 

2. Smooth versions of Oppenheimer-Snyder collapse 

We now consider a family of spacetimes parametrized by a real number a where the 
initial density profiles are also not homogeneous but uniformly converge to a step function 
as 0" — > cxD. These may then be viewed as interpolating between the collapse of a highly 
non-homogeneous dust ball like in the previous example, and Oppenheimer-Snyder collapse. 
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As one would expect, in these examples we will encounter marginally trapped tubes that 
are partially spacelike and partially timelike. 

In particular, consider initial densities of the form^ 

moF(cr) 
Po = ^ 

where F{a) is a complicated function^ chosen so that = limr^oo '^('") is the total dust 
mass. To is the location on the "step" where —dp/dr is a maximum, a characterizes the 
steepness of that step, and erf(a;) is the usual error function. 

In the examples of Figure 4 we take mo as our length scale, choose Tq = 2mo, and consider 
a variety of values of cr. Then graph a) shows how the initial densities converge towards a 
step function with increasing a. In the meantime graph b) shows that for smaller values 
of cr the expansion parameter C is always positive, while for the larger values it starts out 
negative with C ~ —6 while the density is approximately constant. This corresponds to the 
expected value of C in the corresponding OS spacetime. It then diverges to — oo, switches 
to being positive around r = 2mo (when the main part of the step has fallen through), and 
then asymptotes to as the density of dust falling through the horizon also goes to zero. 
As expected, the switch in sign corresponds to the corresponding switch of LnQ{t) — that is 
when the MTT becomes instantaneously tangent to n". 

In c) we also see that in all dynamical horizon asymptotes to R = 2mo as the last 

bits of dust fall through it. However, the MTT behaviour before that time varies greatly. 
For small values of a it continuously increases in area. For large values, we see that there 
are both increasing and decreasing regions. Somewhat confusingly it appears that for a = 2 
there are both increasing and decreasing regions as well, even though b) shows that C > 
everywhere. We will return to examples such as this in section III El, but for now we just 
keep in mind that spacelike surfaces can intersect in non-trivial ways. Here the dynamical 
horizon is demonstrating this as it intersects some of the r = constant surfaces twice. The 
possibility of such foliation effects was also noted in [18]. 

Note that the timelike membranes in this example all go to zero areal radius around 
r = irrrio. This is not surprising as for large values of a, ~ vrmo (eq. (3.18)). Thus, the 
membranes vanish as the first dust shells also collapse to zero areal radius — that is they 
disappear into the density singularity. 

Finally, note that the spacetime diagrams for these spacetimes would be very similar to 
that shown in figure 3. The only significant difference would be that for the spacetimes with 
timelike membranes, the MTT would emerge from/vanish into the singularity as a timelike 
rather than spacelike surface - the slope of the MTT would be greater than 45°. 



erf 



a 



^-1 



(3.26) 



^ Similar initial density profiles can more easily be constructed by means of tanh functions. However, the 
associated mass functions m(r) are awkward expressions in terms of poly logarithms, which are multi- 
branched and need to be treated with care. By defining initial density distributions in terms of error 
functions we avoid such complications. 

3a' 



For those who are interested: F{a) 



27rcr(2cr2 + 3)(1 + erf cr) + 40Fe-'^' (1 + cr^). 
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c) T vs R for MTT evolution 

FIG. 5: Small dust shell 



D. Accretion onto a pre-existing hole 

The next two examples study the accretion of dust shells onto a pre-existing black hole. 
In the first example a black hole will substantially increase its mass while the MTT remains 
null or spacelike everywhere through the evolution. In the second example we study a very 
large matter shell falling into a black hole; here we will again see timelike membranes. 



i. Small dust shell falls into black hole 
We begin with a dust shell of the form 



where measures the total mass of the shell (ie. the asymptotic behaviour of the associated 
mass function), Vo characterizes its "thickness", and a gives its initial position in terms of 
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r = spacelike singularity 




T = 1] = (reflect through this surface) 



FIG. 6: Partial Penrose- Carter diagram showing the accretion of a small dust shell onto a pre- 
existing black hole. To get the full diagram this one should be : 1) continued on its left side with 
the Schwarzschild solution corresponding the the initial mass and 2) reflected through the instant 
of time symmetry. 

To- Again this distribution is Gaussian, but this time it is a shell with peak density at avo. 
We wish to study the accretion of such a shell onto a pre-existing black hole and so excise a 
small region in the interior of the shell spacetime and replace it by a Schwarzschild geometry 
with mass parameter M, putting the junction at some r = f > 2M. From the discussion in 
subsection III A, for r > f the mass function must then be m(r) = M + An po{f)f'^df. 

An example of such a spacetime is shown in Figure 5, where we have chosen f = 2.5 M, 
rrio = M/2, a = 10, and = 10 M. Note that the parameters cannot be chosen arbitrarily in 
this case as shell crossings can easily develop. In fact, a small increase of the mass parameter 
so that nio ~ 0.7 M will be sufficient to cause these. That said, for the parameters that we 
have chosen this does not occur, and we see that the horizon is initially quiescent and begins 
to expand as the dust falls through it. The expansion is spacelike (from figurelllb), C > 0), 
peaks as the largest density of matter crosses the horizon, and then tails off along with the 
infalling dust. Asymptotically the horizon becomes null again with areal radius R = 3M, 
so that the mass function tends to m = |M.^ 

Figure 6 shows the corresponding spacetime diagram. From the spacetime surgery, inside 
the dust shell the spacetime is Schwarzschild with mass M. However, as the dust reaches 
and crosses the MTT, it begins to expand in a spacelike manner and continues to do so 
as long as the dust continues to fall in. Ultimately as that density goes to zero the MTT 
asymptotes to the null event horizon. As in figure 3 the MTT clearly reacts to physical 
events while the event horizon, being teleo logically defined, does not. Note too that as in 
that earlier diagram, the dust distribution is again, for simplicity, shown with an edge. 



^ The asymptotic values will in fact be slightly smaller than that; due to the excision at ?■ = f, the parameter 
TUo is not exactly equal to the mass of the shell. However, in this example the difference is negligible. 
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c) T vs R for MTT evolution 
FIG. 7: Large shell of constant density 



2. Large shell of approximately constant density falls into black hole 

In the previous example, the expansion was everywhere spacelike. It described a fairly 
dramatic situation, a significant expansion of a pre-existing black hole, however there was no 
sign of timelike membranes and the MTTs were isolated or dynamical horizons. We will now 
construct an example with timelike membranes. This will essentially be a smooth version 
of one of the more complicated Oppenheimer-Snyder examples presented in [23]. Thus, we 
build a spacetime in which a large amount of (initially) constant density dust falls into a 
black hole. The initial density is given by 

3m4erf(^)-erf(rg.)) 
^° 47r(r2-ri)(2r2 + 2rir2 + 2r2 + 3M2) ' ^' ^ 

where ri and r2 mark the (approximate) start and end of the shell, rrio is the total mass of 
the shell as read from the mass function. Excision is performed and the interior is replaced 
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by a black hole with mass M, which joins onto the dust exterior at some r = r > 2M. The 
full mass function for r > f will then be m(r) = M + An f! po{f)f'^df. 

For this example we choose ri = 100 M, r2 = 2000 M, and ruo = 600 M. Thus, the mass 
of the hole will increase dramatically during the evolution from m = M to m = 601 M. This 
evolution is shown in figure 7. Note that despite the appearance of the density function, 
it is actually smooth for r > f since the error functions themselves are smooth. Further, 
the density function will "spread" as it falls towards the hole and so the dust will not be of 
constant density as it crosses the horizon. 

That said, we see that the horizon begins to expand in the usual way as the initial 
matter falls into it. Then however, something different happens. Considering evolution with 
respect to r we see that around r = 3000 M , a new marginally trapped surface appears 
at -R = 1202 M and bifurcates into a dynamical horizon and a timelike membrane. The 
dynamical horizon quickly asymptotes to an isolated horizon while the timelike membrane 
contracts and eventually annihilates with the dynamical horizon that grows out to meet it. 

Alternatively if we consider evolution with respect to the initial areal radius r, the MTT 
starts off almost-isolated, becomes dynamical as the first mass falls in, then goes null (with 
— as a tangent) and becomes a timelike membrane that continues to expand as it travels 
backwards in time. Finally, as the last dust falls through, it again goes null (with — n° as a 
tangent), dynamical, and then asymptotes back towards isolation. 

The spacetime diagram for this evolution would be very similar to that shown in figure 
6. The only difference would be that, during its active phase, the MTT would sometimes 
be timelike and so have a slope greater than 45°. 

E. More complicated collapse 

The previous examples display the basic behaviours of marginally trapped tubes — space- 
like expansions, creation from singularities, and timelike contractions/backwards-in-time- 
expansions (depending how one views the evolution). In this section, to get a feel for the 
possible range of evolutions, we will consider examples generated from more complicated 
initial conditions that combine several of these behaviours. 

1. Spacelike expansion/multiple horizons 

This example elaborates on a behaviour that we noted in section III C 2. Namely it shows 
that the appearance of multiple horizons in a given leaf of a spacetime foliation does not 
always signal the existence of timelike membrane sections of the MTT. That is, a dynamical 
horizon as a spacelike surface can interweave with a foliation (of other spacelike surfaces) in 
highly non-trivial ways. 

We will consider matter distributions of the form 

a/i 2 f ^ \ 
I'K'^ror \ To J 

where is an arbitrary reference length scale. When this is integrated to a mass function, 
we see that for any positive integer N there is an amount of mass /i between r = Nirvo/a 
and r = [N + l)nro/a. Thus, physically this corresponds to a series of shells each with the 
same mass (though decreasing density). Again one has to choose parameters with care to 
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c) T vs R for MTT evolution 



FIG. 8: Multiple dust shells fall into a black hole 

avoid shell crossings and/or initial black holes. A particular choice that meets these criteria 
is n = (87r/5)ro and a = 1/4. This is the distribution whose initial configuration and 
evolution is shown in figure 8. 

Then, from Figure 8c), it is clear that r = constant surfaces in this spacetime will contain 
either one or three marginally trapped surfaces (except at turning points of the MTT, where 
they contain two). Further, these will expand and contract in apparently the same kinds of 
ways that we have seen previous example which include timelike membranes. However, an 
examination of Figure 8b), shows that despite this behaviour, the expansion parameter C 
is always greater than zero and so the MTT is everywhere either spacelike or null/isolated. 
Then the apparent contractions/expansions arise simply because the MTT intersects the 
foliation in non-trivial ways. 

In Figure 8c), the horizon goes vertical/null at intervals of i? = 2/i ^ 10 Tq. Each of these 
corresponds to the density going to zero between each shell of mass /i = (87r/5) To ~ 5 and 
so the horizon becoming instantaneously isolated. Finally, note that with a careful choice 
of the parameters, the newly created MTT can absorb shells of this type for arbitrarily long 
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c) T vs R for MTT evolution 



FIG. 9: A more complicated dust collapse. 



periods of time, always alternating between being dynamical and isolated. Thus, a horizon 
can absorb an arbitrarily large amount of mass without ever going timelike. It is the rate 
of absorption, not the total amount that is significant in this regard. 



2. Multiple timelike membranes 



Our final example will demonstrate a spacetime that contains a MTT made up of multiple 
dynamical horizon and timelike membrane regions. The dust density will take the form: 



Po 



a 







f3 + 2cos2 fs^ 



< r < TTTr. 



r > ■KTo 



(3.30) 



where a is a dimensionless constant. The exterior of vrro was excised to avoid negative 
density dust and so violations of the energy conditions. Thus, outside of that radius, the 
geometry will be Schwarzschild with mass M = m{7iro)- 
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Taking a = 1/120 there are neither shell-crossings nor initial black holes in this spacetime. 
The evolution is then shown in the graphs of figure 9. The initial density is irregular and 
gives rise to an MTT which, if we think of evolution as parameterized by r, starts out as a 
dynamical horizon and then alternates back and forth between being timelike and spacelike 
with (non-isolated) null cross-sections separating these regions. From the point of view of 
evolution with respect to r, spacelike slices have anywhere between one and five marginally 
trapped surfaces in this example. However, all of these intricacies are contained within the 
outermost isolated horizon (and here it really is isolated since we cut the density distribution 
at r = nro) and would not be visible to outside observers. 

IV. SCALAR FIELDS 

In the previous section we have seen that the potential MTT behaviours suggested by 
(2.7) and (2.8) are all achieved by the dust spacetimes. Thus, depending on the initial con- 
figuration, the MTT can be any of spacelike, null, or timelike. In this section we will attempt 
the same demonstration for the scalar fields. Thus, we will consider initial configurations of 
scalar fields, evolve them in time, and examine the behaviour of any MTT that forms with 
the help of the expansion parameter C. 

Given that scalar fields are significantly more complicated than dust we will necessarily 
take a numerical rather than an analytic approach. Section IV A will introduce the numer- 
ical model and the subsequent sections will present the results for two different scalar fields 
configurations. These will be analogous to configurations considered in the last section. The 
first, in section IV B, examines the evolution of a initially smooth step-like configuration 
that collapses to form a black hole. The second, in section IV C, studies a shell of scalar 
field which falls into an existing black hole. Due to numerical complexities, these results will 
be less complete than those considered in the last section, but will still demonstrate some 
interesting and complementary behaviours. 

A. Numerical approach 

We will be interested in the scalar fields briefiy considered in section II B 3. The equations 
governing spacetimes containing these fields are generated by the Lagrangian: 



This class of models includes the massive Klein-Gordon field with V = mQcfP' /2 and in the 
following we will restrict ourselves to this case. 

To solve the resulting coupled Einstein-Klein-Gordon equations we work with the stan- 
dard 3+1 approach based on the ADM equations [33] and in particular adopt the techniques 
introduced in [34] to perform the evolutions. Restricting to spherical symmetry, we can then 



There are a few exact scalar field solutions in the literature, see for example [32] and the references listed 
therein. However, they are much more restricted than the Tolman-Bondi solutions (for instance, we are 
not aware of any that are asymptotically flat) and cannot produce the same variety of examples. 




(4.1) 
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study the general evolutions by considering the dynamics of spacetimes with metrics of the 
form 

ds^ = -a^{r, t)dt^ + A{t, r)dr^ + r^B{r, t) {d9^ + sin^ Od'f'^) . (4.2) 

The spherical symmetry implies that all dynamical functions depend only on r and t. a 
is the usual lapse function and for simplicity we impose a vanishing-shift gauge condition; 
equivalently the "time-evolution" vector, d/dt is everywhere orthogonal to the t = constant 
"instantaneous" three-surfaces. 

Focusing on this natural foliation with respect to t, we note that these hypersurfaces have 
intrinsic three-metric and extrinsic curvature : 

hab = A[dr]a[dr]b + Br^riab and (4.3) 

respectively, where flab was defined following equation (3.4). These, together with the value 
of the scalar field, are the variables whose evolution we will study. 

Then, in the usual way we can rewrite the Einstein-Klein- Gordon equations in terms of the 
Hamiltonian and momentum constraints which restrict initial values of these quantities, in 
addition to evolution equations (which preserve the constraints). The actual implementation 
consists of picking an initial configuration for the scalar field, solving the constraints for hab 
and Kab and finally using the evolution equation to obtain the dynamics of the corresponding 
spacetime. For brevity we will not include all of the details here; the interested reader is 
directed to [34] for a more detailed discussion. 

Once we have a spacetime evolution, the next step is to search for an MTT. In this 
case it is natural to search for apparent horizons/marginally trapped surfaces within the 
t = constant slices and this is what we will do. Since we are in spherical symmetry we need 
only consider the null expansions of the r = constant two-surfaces. Then, with 

Ua = -a[dt]a and Sa = \/A[dr]a, (4.5) 

as the forward-in-time and towards-infinity pointing unit normals to the two-surfaces, we 
define ^ 

L = Ua + Sa and = -(«„- s^) (4.6) 

as the outward and inward forward-in-time pointing null normals to those same surfaces. 
As in the dust examples, the exact "normalization" is irrelevant as we are only interested in 
the signs of the ensuing quantities, not their magnitudes. Thus, we can make this specific, 
convenient, choice. Then, it is straightforward to see that: 

^«"> - - ^ + ^•-°^') - H - iiv^ + ^) > (^■«) 

where "Da is the covariant derivative compatible with hab and K = h"'''Kab (as a quick check 
that these expressions are reasonable, note that for Minkowski space where A = B = 1, 
9(^1) = 2/r > and 6'(„) = — 1/r < 0). In the code we evaluate these expressions over 
the whole numerical grid and look for places where 6'(f) changes signs and < 0. The 
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FIG. 10: Scalar field "step" 



corresponding two-spheres are then identified as the marginally trapped two-spheres which 
foliate the MTTs, and their area is 

Aah = ^-ktIhBah (4.9) 

where B^u is the metric function B evaluated at r = tahi the coordinate radius of these 
two-spheres. 

Finally, the expansion parameter C for this surface can be calculated using (2.13). 



B. Scalar field "step" collapses to form a black hole 



Our first example will be analogous to the dust ball collapses considered in section IIIC. 
As for those dust examples, our initial slice will be a moment of time symmetry {Kah = 0) 
and further we will choose our coordinates so that the coordinate radius r will also be the 
areal radius R on that slice (thus B = 1 initially). Then, we specify a step-like configuration 
of the scalar field as shown in figure 10a). 
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Solving the constraint equations to find the initial form of A, we can then integrate the 
data as discussed above to find how the geometry and fields evolve in time. The results are 
shown in figures 10b) and 10c). The first thing to note is that throughout this collapse C > 
and so the MTT is a dynamical horizon. Secondly, as in the earlier examples it asymptotes 
to null as the spacetime settles down to become Schwarzschild. Note from figure 10c) that 
initially the hypersurfaces do not contain an apparent horizon. However, during the collapse 
an apparent horizon appears and then keeps growing until it reaches a constant size. 

C. Scalar field "shell" accretes onto existing black hole 

Our second scalar field example is analogous to the dust shell examples of section HID. 
Thus we will consider an initial "shell" of scalar field that accretes onto an existing black 
hole of mass M. After surgically inserting the black hole we again start out on a slice of 
time symmetry, though in this case for technical reasons do not start with r as the areal 
radius. Instead, on the initial slice B = 1, where S = (l + |^) B. 

Then, as our initial scalar field conditions we consider a 0(0, r) of the form shown in figure 
11a). The corresponding initial intrinsic geometry (essentially only A remains unknown) is 
then found by solving the constraint equations. The results for the corresponding evolution 
are then shown in figures lib) and 11c). On those graphs we again see that C remains 
everywhere non-negative and asymptotes to zero at late times. Note too the apparent 
horizon on the initial hypersurface which grows as the scalar field falls into the hole. 

D. Outlook for scalar fields 

Neither of the preceding examples included timelike membrane regions of the MTT. We 
believe that this is a reflection of the examples that our code has been able to integrate 
rather than a fundamental result. From (2.13) it is clear that a timelike membrane will 
only appear in situations where there is a sufficiently large concentration of the scalar field 
relative to the (inverse) area of the horizon. Our code had difficulty evolving such examples 
and so the lack of timelike membranes is not especially surprising. Physically, one would 
also expect it to be more difficult to obtain such examples for a massive scalar field which, 
in contrast to pressureless dust, resists compression. We expect that future investigations 
will find the appropriate combination of initial conditions needed to generate examples of 
timelike membranes in scalar field spacetimes. For now though, we simply note that our 
examples show timelike evolutions are not the rule for scalar field spacetimes. 

V. CONCLUSIONS AND SPECULATIONS 

In this paper we have seen that dynamical horizons and timelike membranes characterize 
two possible modes of black hole expansion. In the first a black hole/MTT smoothly expands 
as matter falls into it. By contrast, the second occurs when matter densities are high enough 
to force the formation of a new horizon of non-zero area that encloses any already existing 
MTTs. 

Further insight into these two possibilities can be gained by considering the magnitudes 
of the quantities involved. As seen in the discussion of II, assuming that the null energy 
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FIG. 11: Scalar field accretion 



condition holds, the signature of an MTT is determined by the relative magnitudes of 1/^1 
and 2Tafe^"n^, or, specializing to (pressureless) dust, 1/A and p. Converting into physical 
units, it is straightforward to see that for a Schwarzschild black hole of mass M = hMq, the 
inverse horizon area corresponds to a density 



Pa 



(5.1) 



where Mq = 2 x 10^^ g is the approximate mass of the Sun and c and G are respectively 
the speed of light and the gravitational constant. Then, for a solar mass black hole, pa ~ 
10^^ g/cm^ which is about an order of magnitude higher than the density of a neutron star. 
By contrast, for a supermassive black hole of mass 10^ Mq, one has p^ ^ 1 g/cm^, the 
density of water. 

Though strictly speaking these results only apply in situations of spherical symmetry, it 
seems fairly safe to use them to draw some more general conclusions. Specifically, it is likely 
that both modes of expansion are not only mathematically possible as we have seen in this 
paper, but also both occur in physical situations. For small black holes, even if the lack 
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of spherical symmetry changes these estimates by several orders of magnitude, it appears 
that dynamical horizon spacelike expansions are probably the dominant mode in all but 
the most extreme situations, such as black hole and/or neutron star collisions. Numerical 
studies support this contention as in such extreme situations the occurrence of multiple 
horizons appears to be generic [21, 35] (though also relatively unstudied in a systematic 
way, since in most studies it is the exterior spacetime that is of interest and so the interior of 
the outermost apparent horizon is excised and thrown away). By contrast for supermassive 
black holes, it is likely that horizon jumps/TLMs are much more common. 

These examples also suggest other (possible) properties of MTTs. First, dynamical hori- 
zons appear to only originate either out of singularities (the density singularities of 111 C) or 
as part of a dynamical horizon/TLM pair (as in HIE 2). Equivalently, in all of our examples 
there is just one MTT associated with each black hole. This originates in a singularity and 
then weaves backwards and forwards in time, always expanding in area (relative to a folia- 
tion parameter that monotonically increases as one moves away from the singularity). This 
suggests that a similar result may be true away from spherical symmetry — the multiple 
horizons/jumps seen in numerical studies of black hole collisions may actually all be part of 
a single MTT that weaves backwards and forwards in time. 

In our examples it is also true that TLMs and dynamical horizons always occur in such a 
way that a causal signal originating from the MTT would never be detectable by sufficiently 
far-away observers. This is consistent with a gravitational confinement theorem due to 
Israel [29, 36] which states that if the weak energy condition holds (as is the case in all our 
examples), a trapped two-sphere can be extended to a spacelike three-cylinder foliated by 
trapped two-spheres of constant area. Assuming reasonably regular spacetimes, this three- 
cylinder will act as a permanent one-way membrane for causal effects. Even though an 
observer "inside" an MTT can escape that trapped region, he will not be able to send signals 
beyond the areal radius of the two-sphere on which the MTT was first crossed. We note 
that Israel's confinement theorem is quite general; it does not assume spherical symmetry or 
even asymptotic fiatness. Moreover, in appropriate asymptotically fiat spacetimes, trapped 
surfaces must necessarily be contained within event horizons and so unable to send causal 
signals to null infinity [1, 2]. Our asymptotically fiat examples are also consistent with the 
latter: the outermost parts of MTTs are dynamical horizons which asymptote to a null 
surface with some finite area. Consequently, in physically realistic spacetimes, any TLMs 
that may be present will be hidden from observers far from a black hole, also in the absence 
of spherical symmetry. 

In summary, even though the examples and (dust) calculations that we have seen in 
this paper are quite simple we believe that they are extremely useful in forming a correct 
intuition about the behaviour of MTTs during general black hole evolution. In particular 
they provide a convenient testing ground for ideas about these evolutions which is much 
simpler than full numerical simulations of black hole collisions and yet still significantly 
richer (and more realistic) than the heretofore studied analytical examples (Vaidya and 
Oppenheimer-Snyder). As such they should be useful in, among other things, suggesting 
possible extensions of the recent mathematical investigations [17, 18] of MTT properties. 
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